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Abstract 

The fc-ellipse is the plane algebraic curve consisting of all points whose sum 
of distances from k given points is a fixed number. The polynomial equation 
defining the £>ellipse has degree 2 k if k is odd and degree 2 k — ( fc , 2 ) if k is even. 
We express this polynomial equation as the determinant of a symmetric matrix 
of linear polynomials. Our representation extends to weighted /c-ellipses and k- 
ellipsoids in arbitrary dimensions, and it leads to new geometric applications of 
semidefinite programming. 



1 Introduction 



y - vi 

d — x + Ui 



The circle is the plane curve consisting of all points (x, y) whose distance from a given 
point (ui,vx) is a fixed number d. It is the zero set of the quadratic polynomial 

/ \ i , + x — ui 

Pi(x,y) = det 

[ y-vi 

The ellipse is the plane curve consisting of all points (x, y) whose sum of distances from 
two given points (ui, v{) and (u 2 , v 2 ) is a fixed number d. It is the zero set of 

d + 2x — u\ — ui y — Vi y — v% 

y — v\ d + u\ — U2 V — V2 

y — v 2 d — u\ + U2 y — v± 

y — f 2 y — vi d — 2x + u\ + u 2 

In this paper we generalize these determinantal formulas for the circle and the ellipse. 
Fix a positive real number d and fix k distinct points (ui, v i), (112, v 2 ), ■ ■ ■ , (v,k, Vk) in 
IR 2 . The k-ellipse with foci (ui,Vi) and radius d is the following curve in the plane: 

k 

(x, y ) em 2 y/(z - ^) 2 + (y- = d } ■ ('•>) 



P2(x,y) 



det 



(2) 
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Figure 1: A 3-ellipse, a 4-ellipse, and a 5-ellipse, each with its foci. 



The fc-ellipse is the boundary of a convex set £ k in the plane, namely, the set of 
points whose sum of distances to the k given points is at most d. These convex sets are 
of interest in computational geometry [Tra06] and in optimization, e.g. for the Fermat- 
Weber facility location problem |Baj88| ICT901 IKul77} ISek991 IWel37] . In the classical 



literature (e.g. [Stu84] ) . fc-ellipses are known as Tschirnhaus'sche Eikurven |N agy50| . 
Indeed, they look like "egg curves" and they were introduced by Tschirnhaus in 1686. 

We are interested in the irreducible polynomial p k (x,y) that vanishes on the k- 
ellipse. This is the unique (up to sign) polynomial with integer coefficients in the 
unknowns x and y and the parameters d, U\,V\, . . . ,u k ,v k . By the degree of the k- 
ellipse we mean the degree of Pk(x,y) in x and y. To compute it, we must eliminate 
the square roots in the representation ([3]). Our solution to this problem is as follows: 



Theorem 1.1. The k-ellipse has degree 2 k if k is odd and degree 2 k — ( fc , 2 ) if k 



is even. 



Its defining polynomial has a determinantal representation 

Pk(x,y) = det(x • A k + y- B k + C k ) (4) 

where A k , B k , C k are symmetric 2 k x 2 k matrices. The entries of A k and B k are integer 
numbers, and the entries of C k are linear forms in the parameters d, U\,V\, . . . , u k , v k . 

For the circle (k = 1) and the ellipse (k = 2), the representation (T4j) is given by the 
formulas (CEJ) and (T5]). The polynomial p^{x,y) for the 3-ellipse is the determinant of 

d-\-3x — u \ — it 2 —u^y — vi y — V2 y — v-$ 

y — v\ d-\-x-\-u\ —U2 — li 3 y — V2 y — 

y — V2 d-\-x — ui -\-U2 —u^y — vi y — U3 

y — V2 y~ v l d — x-\-ui +112 — 1A3O y — l>3 

y — -U3 d-\-x — u\ — 1A2+1A3 y — v\ y — 1>2 

y— U3 y— v\ d—x-\-u\ —112+^3 y—^i 

y — U3 y — V2 d — x — ui-\-U2 J rU3 y~ v l 

y — U3 y—V2 y~ v i rf— 3a; +u, 1+^2+1^3. 

The full expansion of this 8 x 8-determinant has 2, 355 terms. Next, the 4-ellipse is a 
curve of degree ten which is represented by a symmetric 16 x 16- matrix, etc.... 

This paper is organized as follows. The proof of Theorem 1 1.1 1 will be given in Section 
2. Section 3 is devoted to geometric aspects and connections to semidefmite program- 
ming. While the fc-ellipse itself is a convex curve, its Zariski closure { p k (x,y) = 0} 
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Figure 2: The Zariski closure of the 5-ellipse is an algebraic curve of degree 32. 



has many extra branches outside the convex set They are arranged in a beautiful 
pattern known as a Helton- Vinnikov curve [H V03j . This pattern is shown in Figure [2] 
for k = 5 points. In Section 4 we generalize our results to higher dimensions and to the 
weighted case, and we discuss the computation of the Fermat- Weber point of the given 
points {ui, Vi). A list of open problems and future directions is presented in Section 5. 

2 Derivation of the Matrix Representation 

We begin with a discussion of the degree of the fc-ellipse. 

Lemma 2.1. The defining polynomial of the k-ellipse has degree at most 2 k in the 
variables (x,y) and it is monic of degree 2 k in the radius parameter d. 

Proof. We claim that the defining polynomial of the fc-ellipse can be written as follows: 
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Obviously, the right hand side vanishes on the fc-ellipse. The following Galois theory 
argument shows that this expression is a polynomial and that it is irreducible. Consider 
the polynomial ring R = Q[x, y, d, u\, vi, . . . , u k , v k ]. The field of fractions of R is the 
field K = Q(x, y,d,Ui,Vi, . . . ,Uk, ujfe) of rational functions in all unknowns. Adjoining 
the square roots in ([3]) to K gives an algebraic field extension L of degree 2 k over K. 
The Galois group of the extension L/K is (Z/2Z) fc , and the product in (j3J) is over 
the orbit of the element d — Yli=i V ( x ~ Ui ) 2 + (y ~ Vi ) 2 °^ ^ under the action of the 
Galois group. Thus this product in (jSJ) lies in the ground field K. Moreover, each factor 
in the product is integral over R, and therefore the product lies in the polynomial ring 
R. To see that this polynomial is irreducible, it suffices to observe that no proper 
subproduct of the right hand side in lies in the ground field K. □ 

The statement degree at most 2 k is the crux in Lemma 12.11 Indeed, the degree in 
(x, y) and can be strictly smaller than 2 k as the case of the classical ellipse (k = 2) 
demonstrates. When evaluating the product ((SJ) some unexpected cancellations may 
occur. This phenomenon happens for all even k, as we shall see later in this section. 

We now turn to the matrix representation promised by Theorem 11.11 We recall 
the following standard definition from matrix theory (e.g., [HJ94j ). Let A be a real 
m x m-matrix and B a real n x n-matrix. The tensor sum of A and B is the mn x mn 
matrix A® B := A® I n + I m ® B . The tensor sum of square matrices is an associative 
operation which is not commutative. For instance, for three matrices A, B, C we have 



Here (g> denotes the tensor product, which is also associative but not commutative. 
Tensor products and tensor sums of matrices are also known as Kronecker products and 
Kronecker sums [Bel97, HJ94J. Tensor sums of symmetric matrices can be diagonalized 
by treating the summands separately: 

Lemma 2.2. Let M 1; . . . , be symmetric matrices, Ui, . . . , orthogonal matrices, 
and Ai, . . . , Afc diagonal matrices such that Mi = Ui • Aj • Uj for i = 1, . . . , k. Then 



In particular, the eigenvalues of the tensor sum M\ © M2 © • • • © are the sums 
Ai + A2 + • • • + Afc where Ai is any eigenvalue of Mi, A2 is any eigenvalue of M2, etc. 

The proof of this lemma is an exercise in (multi)-linear algebra. We are now pre- 
pared to state our formula for the explicit determinantal representation of the fc-ellipse. 

Theorem 2.3. Given points (ui,Vi), . . . , (uk,Vk) in WL 2 , we define the 2 k x 2 k matrix 



A®B®C 



A® I® I + I®B®I + I®I®C. 



(C/i <g> • • • <g> U k ) T ■ (Mi ® • • • © M k ) ■ (C/i <g> • • • ® U k ) = Ai © ■ • • © Afc. 



L k {x,y) 



d ■ I 2 k + 



x — u\ y — v\ 
y — v\ —x + u\ 



© 



© 



x - u k y - v k 
y-v k -x + u k 



(6) 
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which is affine in x,y and d. Then the k-ellipse has the determinantal representation 

Pk(x,y) = det L k (x,y). (7) 
The convex region bounded by the k-ellipse is defined by the following matrix inequality: 

£k = {(x,y)eR 2 : L k (x,y) h 0}. (8) 
Proof. Consider the 2x2 matrix that appears as a tensor summand in (jSJ): 



y-Vi 



y-v t 

-X + Ui 



A computation shows that this matrix is orthogonally similar to 



[x 



U; 



+ (y 



o 



Ui 



+ (y 



These computations take place in the field L which was considered in the proof of 
Lemma 12.11 above. Lemma 12.21 is valid over any field, and it implies that the matrix 
Lk(x,y) is orthogonally similar to a 2 k x 2 k diagonal matrix with diagonal entries 



d + ^(-i)-V(x- Ui ) 2 + (y-^) 2 , <n e {o, i}. 



(9) 



i=l 



The desired identity ([7]) now follows directly from (j3J) and the fact that the determinant 
of a matrix is the product of its eigenvalues. For the characterization of the convex 
set £k, notice that positive semidefiniteness of Lk(x,y) is equivalent to nonnegativity 
of all the eigenvalues Qj. It suffices to consider the smallest eigenvalue, which equals 

k 

d - ^2 V( x ~ u d 2 + (y- v i) 2 - 



Indeed, this quantity is nonnegative if and only if the point (x, y) lies in 



□ 



Proof of Theorem \l.l[ The assertions in the second and third sentence have just been 
proved in Theorem 12.31 What remains to be shown is the first assertion concerning the 
degree of pk(x,y) as a polynomial in (x,y). To this end, we consider the univariate 
polynomial g(t) := pk(t cos 6 , t sin 8) where 9 is a generic angle. We must prove that 



deg t (<?(*)) 



(fc/2) 



if k is odd, 
if k is even. 



The polynomial g(t) is the determinant of the symmetric 2 k x 2 fc -matrix 



Lk(t cos 9, t sin 1 



cos# 
sin 9 



sin# 

cos 9 



cos 9 
sin 9 



sin 9 

cos 9 



Ch. 



(10) 
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The matrix does not depend on t. We now define the 2 k x 2 k orthogonal matrix 



U := V< 

k times 

and we note the matrix identity 



where V 



V 1 ■ 



cos 9 
sin 9 



sin 9 

cos 9 



■ V 



cos(0/2) 
sin (0/2) 



1 
-1 



•sin(0/2)' 
cos(0/2) 



Pre- and post-multiplying (|T0|) by £/ T and U, we find that 
f/ T • L fc (tcos0,tsin0) ■ U = t 



1 
-1 





"1 " 






-1 


) 



U T ■ c, ■ u. 



E k 

The determinant of this matrix is our univariate polynomial git)- The matrix is a 
diagonal matrix of dimension 2 k x 2 k . Its diagonal entries are obtained by summing k 
copies of —1 or +1 in all 2 k possible ways. None of these sums are zero when k is odd, 
and precisely ( fe , 2 ) of these sums are zero when k is even. This shows that the rank 

of Ek is 2 k when k is odd, and it is 2 k — (/j 2 ) when is even. We conclude that the 

univariate polynomial g{t) = det ( t ■ Ek + U T CkU) has the desired degree. □ 



3 More Pictures and Some Semidefinite Aspects 

In this section we examine the geometry of the fc-ellipse, we look at some pictures, and 
we discuss aspects relevant to the theory of semidefinite programming. In Figure Q] 
several fc-ellipses are shown, for k = 3, 4, 5. One immediately observes that, in contrast 
to the classical circle and ellipse, a fc-ellipse does not necessarily contain the foci in its 
interior. The interior Ek of the /c-ellipse is a sublevel set of the convex function 

k 

(x, y ) h- J2V^-^) 2 + (y-^) 2 - (ii) 

This function is strictly convex in any region excluding the foci {(iii, v «)}f =1 , provided 
the foci are not collinear [Sek99j. This explains why the /c-ellipse is a convex curve. In 
order for Sk to be nonempty, it is necessary and sufficient that the radius d be greater 
than or equal to the global minimum d+ of the convex function ffTTl) . 

The point (x±, y±) at which the global minimum d* is achieved is called the Fermat- 
Weber point of the foci. This point minimizes the sum of the distances to the k 
given points (wj,t>j), and it is of importance in the facility location problem. See 
|Baj88[ ICT901 [Kul77] . and |Stu84j for a historical reference. For a given set of foci, we 
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Figure 3: A pencil of 3-ellipses with fixed foci (the three black dots) and different radii. 
These convex curves are always smooth unless they contain one of the foci. 

can vary the radius d, and this results in a pencil of confocal fc-ellipses, as in Figure EJ 
The sum of distances function f jTTj) is differentiable everywhere except at the (ui,Vi), 
where the square root function has a singularity. As a consequence, the fc-ellipse is a 
smooth convex curve, except when that curve contains one of the foci. 

An algebraic geometer would argue that there is more to the fc-ellipse than meets 
the eye in Figures [1] and [31 We define the algebraic k-ellipse to be the Zariski closure 
of the fc-ellipse, or, equivalently, the zero set of the polynomial pk{x,y). The algebraic 
fc-ellipse is an algebraic curve, and it can be considered in either the real plane R 2 , in 
the complex plane C 2 , or (even better) in the complex projective plane P^. 

Figure [2] shows an algebraic 5-ellipse. In that picture, the actual 5-ellipse is the 
tiny convex curve in the center. It surrounds only one of the five foci. 

For a less dizzying illustration see Figure HI That picture shows an algebraic 3- 
ellipse. The curve has degree eight, and it is given algebraically by the 8 x 8-determinant 
displayed in the Introduction. We see that the set of real points on the algebraic 3- 
ellipse consists of four ovals, corresponding to the equations 

y/(x- Mi) 2 + (y- v i) 2 ±y/(x- u 2 ) 2 + (y- v 2 ) 2 ±y/(x- u 3 ) 2 + (y- v 3 ) 2 = d. 

Thus Figure H] visualizes the Galois theory argument in the proof of Lemma 12.11 

If we regard the radius d as an unknown, in addition to the two unknowns x and y, 
then the determinant in Theorem 11.11 specifies an irreducible surface {pk(x,y,d) = 0} 
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Figure 4: The Zariski closure of the 3-ellipse is an algebraic curve of degree eight. 



in three-dimensional space. That surface has degree 2 k . For an algebraic geometer, 
this surface would live in complex projective 3-space CP 3 , but we are interested in its 
points in real affine 3-space M 3 . Figure shows this surface for k = 3. The bowl-shaped 
convex branch near the top is the graph of the sum of distances function (ITT)) , while 
each of the other three branches is associated with a different combination of signs in 
the product ([5]). The surface has a total of 2 k = 8 branches, but only the four in the 
half-space d > are shown, as it is symmetric with respect to the plane d = 0. Note 
that the Fermat- Weber point (x*, y*, d*) is a highly singular point of the surface. 

The time has now come for us to explain the adjective "semidefmite" in the title 
of this paper. Semidefinite programming (SDP) is a widely used method in convex op- 
timization. Introductory references include VBiHil IWSV00] . An algebraic perspective 
was recently given in [NRS06] . The problem of SDP is to minimize a linear functional 
over the solution set of a linear matrix inequality (LMI). An example of an LMI is 



Here Ck is the matrix gotten from Ck by setting d — 0, so that Ck = d- I 2 k + Ck- If d is 
a fixed positive real number then the solution set to the LMI (IT2")) is the convex region 
Sk bounded by the fc-ellipse. If d is an unknown then the solution set to the LMI f[T21) 
is the epigraph of fTTTj) . or, geometrically, the unbounded 3-dimensional convex region 
interior to the bowl-shaped surface in Figure [5j The bottom point of that convex region 
is the Fermat- Weber point (x*, y*, <i*), and it can be computed by solving the SDP 



x ■ A k + y ■ B k + d ■ I 2 k + C k y 0. 



(12) 



Minimize d subject to (1T21 . 



(13) 
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Figure 5: The irreducible surface {p^(x, y, d) = 0}. Taking horizontal slices gives the 
pencil of algebraic 3-ellipses for fixed foci and different radii d, as shown in Figure El 



Similarly, for fixed radius d, the fc-ellipse is the set of all solutions to 

Minimize ax + (3y subject to (1T21) 



(14) 



where a, (3 run over R. This explains the term semidefinite representation in our title. 

While the Fermat- Weber SDP (fT3l) has only three unknowns, it has the serious dis- 
advantage that the matrices are exponentially large (size 2 k ). For computing (sc*, y*, d*) 
in practice, it is better to introduce slack variables d\, g?2> • • • , d/~, and to solve 



Minimize di subject to 



di + x-cii y-k 
y-bi di-x + ai 



y (i = l,...,k). (15) 



This system can be written as single LMI by stacking the 2 x 2-matrices to form a block 
matrix of size 2k x 2k. The size of the resulting LMI is linear in k while the size of the 
LMI (fT51) is exponential in k. If we take the LMI (|15p and add the linear constraint 
di -\- d 2 + ■ ■ ■ + dk = d, for some fixed d > 0, then this furnishes a natural and concise 
lifted semidefinite representation of our fc-ellipse. Geometrically, the representation 
expresses Sk as the projection of a convex set defined by linear matrix inequalities in a 
higher- dimensional space. Theorem 11.11 solves the algebraic LMI elimination problem 
corresponding to this projection, but at the expense of an exponential increase in size, 
which is due to the exponential growing of degrees of /c-ellipses. 
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Our last topic in this section is the relationship of the /c-ellipse to the celebrated 
work of Helton and Vinnikov [HV03j on LMI representations of planar convex sets, 
which led to the resolution of the Lax conjecture in |LPR05j . A semialgebraic set 
in the plane is called rigidly convex if its boundary has the property that every line 
passing through its interior intersects the Zariski closure of the boundary only in real 
points. Helton and Vinnikov [HV03t Thm. 2.2] proved that a plane curve of degree 
d has an LMI representation by symmetric d x d matrices if and only if the region 
bounded by this curve is rigidly convex. In arbitrary dimensions, rigid convexity holds 
for every region bounded by a hypersurface that is given by an LMI representation, 
but the strong form of the converse, where the degree of the hypersurface precisely 
matches the matrix size of the LMI, is only valid in two dimensions. 

It follows from the LMI representation in Theorem 11.11 that the region bounded by 
a /c-ellipse is rigidly convex. Rigid convexity can be seen in Figures H] and [2j Every 
line that passes through the interior of the 3-ellipse intersects the algebraic 3-ellipse 
in eight real points, and lines through the 5-ellipse meet its Zariski closure in 32 real 
points. Combining our Theorem 11.11 with the Helton- Vinnikov Theorem, we conclude: 



Corollary 3.1. The k-ellipse is rigidly convex. If k is odd, it can be represented by an 

I&/2) 



LMI of size 2 k , and if k is even, it can be represented as by LMI of size 2 k — (A) . 



We have not found yet an explicit representation of size 2 k — (A) when k is even and 
k > 4. For the classical ellipse (k — 2), the determinantal representation (j2J) presented 
in the Introduction has size 4x4, while Corollary 13.11 guarantees the existence of a 2 x 2 
representation. One such representation of the ellipse with foci (iti,i>i) and (u2,v 2 ) is: 



(d 2 + (ui-U2)(2x-ui-u 2 ) + {vi-v 2 )(2y-vi-V2)) ■ 



h + 2d- 



x - u 2 y - v 2 
y - V2 -x + u 2 



y 0. 



Notice that in this LMI representation of the ellipse, the matrix entries are linear in x 
and y, as required, but they are quadratic in the radius parameter d and the foci ttj, Uj. 
What is the nicest generalization of this representation to the /c-ellipse for k even ? 



4 Generalizations 

The semidefmite representation of the /c-ellipse we have found in Theorem 11.11 can be 
generalized in several directions. Our first generalization corresponds to the inclusion of 
arbitrary nonnegative weights for the distances, while the second one extends the results 
from plane curves to higher dimensions. The resulting geometric shapes are known 
Tschirnhaus'sche Eifldchen (or "egg surfaces") in the classical literature |Nagy50| . 
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4.1 Weighted /c-ellipse 

Consider k points (ui,Vx), . . . , (v,k,Vk) in the real plane M 2 , a positive radius d, and 
positive weights w±, . . . , w k . The weighted k- ellipse is the plane curve defined as 

k 

( i,|/)6l 2 : Wi ■ a/(x - Mj) 2 + (y - v { ) 2 = d 
i=i 

where Wi indicates the relative weight of the distance from (x, y) to the i-th focus 
(ui, Vi). The algebraic weighted k-ellipse is the Zariski closure of this curve. It is the zero 
set of an irreducible polynomial p k (x,y) that can be constructed as in equation ([5]). 
The interior of the weighted fc-ellipse is the bounded convex region 



£ k {w) :-- 



GR 2 : Y^ Wi - ^{x-Uif + iy-v.f < d J 



The matrix construction from the unweighted case in generalizes as follows: 



Ll"(x, y) := d ■ I 2 k + Wi 



x — U\ y — vi 
y — v i —x + ui 



x-u k y -v k 
y-v k -x + u k 



(16) 



Each tensor summand is simply multiplied by the corresponding weight. The following 
representation theorem and degree formula are a direct generalization of Theorem 12.31 

Theorem 4.1. The algebraic weighted k-ellipse has the semidefinite representation 

p w k (x,y) = det Lt(x,y), 

and the convex region in its interior satisfies 

S k (w) = {(x,y) GR 2 : L%(x,y) h 0}. 

The degree of the weighted k-ellipse is given by 

deg Pk »{x,y) = 2 k -\V(w)\, 

where V(w) = {5 G {-1, l} fc : Eti = 0}- 

Proof. The proof is entirely analogous to that of Theorem 12. 3t □ 

A consequence of the characterization above is the following cute complexity result. 

Corollary 4.2. The decision problem "Given a weighted k-ellipse with fixed foci and 
positive integer weights, is its algebraic degree smaller than 2 k ? " is NP-complete. 

Proof. Since the number partitioning problem is NP-complete [G J79] , this follows from 
the degree formula in Theorem 14. 1L 

□ 
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4.2 /c-Ellipsoids 

The definition of a /c-ellipse in the plane can be naturally extended to a higher- 
dimensional space to obtain /c-ellipsoids. Consider k fixed points ui,...,Ufc in M. n , 
with Uj = (uii, Ui2, ■ ■ ■ , Ui n ), the k-ellipsoid in R n with these foci is the hypersurface 



{ x G R n : IK - x ll = d} = \xeW n :J2 



t=l 



\ 



3=1 



d \ . (17) 



This hypersurface encloses the convex region 



or 



{xGR" : ^||ui-x|| < d). 



i=l 



The Zariski closure of the fc-ellipsoid is defined by an irreducible polynomial p£(x) = 
Pk(xi, x 2 , ■ ■ ■ , x n ). By the same reasoning as in Section 2, we can prove the following: 

Theorem 4.3. The defining irreducible polynomial p£(x) o/ t/ie k-ellipsoid is monic 
of degree 2 k in the parameter d, it has degree 2 k in x i/ A; is odd, and it has degree 
2 k - (A) i/ fc is even. 



We shall represent the polynomial p„(x) as a factor of the determinant of a sym- 
metric matrix of affine-linear forms. To construct this semidefinite representation of 
the ^-ellipsoid, we proceed as follows. Fix an integer m > 2. Let U,(x) be any sym- 
metric mxm-matrix of rank 2 whose entries are affine-linear forms in x, and whose 
two non-zero eigenvalues are ±||uj — x||. Forming the tensor sum of these matrices, as 
in the proof of Theorem 12.31 we find that p^. (x) is a factor of 



det (d-I m k + Ui(x)©U 2 (x) 



U fc (x)). 



18) 



However, there are many extraneous factors. They are powers of the irreducible poly- 
nomials that define the /c'-ellipsoids whose foci are subsets of {u 1; u 2 , . . . , u^}. 

There is a standard choice for the matrices U»(x) which is symmetric with respect 
to permutations of the n coordinates. Namely, we can take m = n + 1 and 



Ui(x) 



xi—ua x 2 -u i2 
x\ — Ua 
x 2 —u i2 



U; 











X n Hi 

o 









However, in view of the extraneous factors in ([TBI , it is desirable to replace these 
by matrices of smaller size m, possibly at the expense of having additional nonzero 
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eigenvalues. It is not hard to see that n is the smallest possible matrix size in a 
symmetric determinant representation. To see this, assume d 2 — YTj=i( x j ~ u ij) 2 = 
det(A + xiAi + - ■ -+x n A n ) where A4 are all constant symmetric matrices of size mxm. 
Then it must hold m > n. Otherwise, if m < n, the first row of Aq + x\A\ + • • ■ + x n A n 
would vanish on some affine subspace of R™ and so did d 2 — YTj=i( x j ~ u ij) 2 1 but this 
is not possible since {x G M" : d 2 — Y^j=i(. x j ~~ u ij) 2 = 0} * s compact. However, if 
we allow extraneous factors, then the smallest possible value of m might drop. These 
questions are closely related to find the determinant complexity of a given polynomial, 
as discussed in |MR04] . Note that in the applications to complexity theory considered 
there, the matrices of linear forms need not be symmetric. 

5 Open questions and further research 

The /c-ellipse is an appealing example of an object from algebraic geometry. Its defi- 
nition is elementary and intuitive, and yet it serves well in illustrating the intriguing 
interplay between algebraic concepts and convex optimization, in particular semidefi- 
nite programming. The developments presented in this paper motivate many natural 
questions. For most of these, to the best of our knowledge, we currently lack definite 
answers. Here is a short list of open problems and possible topics of future research. 

Singularities and genus Both the circle and the ellipse are rational curves, i.e., 
have genus zero. What is the genus of the (projective) algebraic /c-ellipse? The first 
values, from k = 1 to k = 4, are 0,0,3,6. What is the formula for the genus in general? 
The genus is related to the class of the curve, i.e. the degree of the dual curve, and this 
number is the algebraic degree |NRS 06j of the problem (fl4"I) . Moreover, is there a nice 
geometric characterization of all (complex) singular points of the algebraic fc-ellipse? 

Algebraic degree of the Fermat- Weber point The Fermat- Weber point 
is the unique solution of an algebraic optimization problem, formulated in (fl3l or ([15]) . 
and hence it has a well-defined algebraic degree over Q(«i, V\, . . . , u^, ffc). However, 
that algebraic degree will depend on the combinatorial configuration of the foci. For 
instance, in the case k = 4 and foci forming a convex quadrilateral, the Fermat- Weber 
point lies in the intersection of the two diagonals |Baj88|, and therefore its algebraic 
degree is equal to one. What are the possible values for this degree? Perhaps a possible 
approach to this question would be to combine the results and formulas in |NRS06] 
with the semidefinite characterizations obtained in this paper. 

Reduced SDP representations of rigidly convex curves A natural question 
motivated by our discussion in Section [3] is how to systematically produce minimal 
determinantal representations for a rigidly convex curve, when a non-minimal one is 
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available. This is likely an easier task than finding a representation directly from the 
defining polynomial, since in this case we have a certificate of its rigid convexity. 
Concretely, given real symmetric n x n matrices A and B such that 

p(x, y) = det(A ■ x + B ■ y + I n ) 

is a polynomial of degree r < n, we want to produce r x r matrices A and B such that 

p(x, y) = det(A ■ x + B ■ y + I r ). 

The existence of such matrices is guaranteed by the results in [HV031 [LPR 05J . In fact, 
explicit formulas in terms of theta functions of a Jacobian variety are presented in 
[HV03J. But isn't there a simpler algebraic construction in this special case? 



Elimination in semidefinite programming The projection of an algebraic variety 
is (up to Zariski closure, and over an algebraically closed field) again an algebraic 
variety. That projection can be computed using elimination theory or Grobner bases. 
The projection of a polyhedron into a lower-dimensional subspace is a polyhedron. 
That projection can be computed using Fourier-Motzkin elimination. In contrast to 
these examples, the class of feasible sets of semidefinite programs is not closed under 
projections. As a simple concrete example, consider the convex set 



[x,y,t) e 



1 x — t 
x — t y 



>- 0, t > 



Its projection onto the (x, ?/)-plane is a convex set that is not rigidly convex, and hence 
cannot be expressed as {(x,y) e M 2 : Ax + By + C >z 0}. In fact, that projection is not 
even basic semialgebraic. In some cases, however, this closure property nevertheless 
does hold. We saw this for the projection that transforms the representation (|T5l) 
of the fc-ellipse to the representation (|T3l) . Are there general conditions that ensure 
the semidefinite representability of the projections? Are there situations where the 
projection does not lead to an exponential blowup in the size of the representation? 



Hypersurfaces defined by eigenvalue sums Our construction of the (weighted) 
fc-ellipsoid as the determinant of a certain tensor sum has the following natural general- 
ization. Let Ui(x), U2(x), . . . , Ufc(x) be any symmetric m x m- matrices whose entries 
are affine-linear forms in x = (xi, X2, ■ ■ ■ , x n ). Then we consider the polynomial 

p(x) = det(Ui(x)©U 2 (x)©---©U fc (x)). (19) 

We also consider the corresponding rigidly convex set 

{x£R" : Ui(x) © U 2 (x) © ■ ■ ■ © Ua(x) h }. 
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The boundary of this convex set is a hypersurface whose Zariski closure is the set of 
zeroes of the polynomial p(x). It would be worthwhile to study the hypersurfaces of 
the special form (lT9~j) from the point of view of computational algebraic geometry. 

These hypersurfaces specified by eigenvalue sums of symmetric matrices of linear 
forms have a natural generalization in terms of resultant sums of hyperbolic polynomi- 
als. For concreteness, let us take k = 2. If p(x) and q(x) are hyperbolic polynomials 
in n unknowns, with respect to a common direction e in R", then the polynomial 

(p©g)(x) := Res t (p(x - te), g(x + te)) 

is also hyperbolic with respect to e. This construction mirrors the operation of taking 
Minkowski sums in the context of convex polyhedra, and we believe that it is funda- 
mental for future studies of hyperbolicity in polynomial optimization [LPR05t IRen06| . 
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